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Dark energy must cluster in order to be consistent with the equivalence principle. The background 
evolution can be effectively modeled by either a scalar field or by a barotropic fluid. The fluid model 
can be used to emulate perturbations in a scalar field model of dark energy, though this model breaks 
down at large scales. In this paper we study evolution of dark energy perturbations in canonical 
scalar field models: the classes of thawing and freezing models. The dark energy equation of state 
evolves differently in these classes. In freezing models, the equation of state deviates from that 
of a cosmological constant at early times. For thawing models, the dark energy equation of state 
remains near that of the cosmological constant at early times and begins to deviate from it only 
at late times. Since the dark energy equation of state evolves differently in these classes, the dark 
energy perturbations too evolve differently. In freezing models, since the equation of state deviates 
from that of a cosmological constant at early times, there is a significant difference in evolution 
of matter perturbations from those in the cosmological constant model. In comparison, matter 
perturbations in thawing models differ from the cosmological constant only at late times. This 
difference provides an additional handle to distinguish between these classes of models and this 
difference should manifest itself in the ISW effect. 

PACS numbers: 



I. INTRODUCTION 



Various observation have confirmed that the expansion 
rate of the universe is accelerating [l|. These observa- 
tions include those of Supernova type la 0] , observations 
of Cosmic Microwave Background [3, 0] and large scale 
structure [5( . The accelerated expansion of the universe 
can be explained by introducing a cosmological constant 
A in the Einstein's equation [y, 0]. However, the cos- 
mological constant model is plagued by the fine tuning 
problem Q. This has motivated the study of dark en- 
ergy models to explain the current accelerated expansion 
of the universe (for reviews see @). An alternative to 
the cosmological constant model is to assume that this 
accelerated expansion is driven by a canonical scalar field 
with a potential V((f), namely the quintessence field 0- 
[T3 |. There exists another class of string theory inspired 
scalar field dark energy models known as tachyon models 
[TH [l9j and there are models which allow w < — 1 are 
known as phantom models [20T |. Phantom type dark en- 
ergy can also be realized in a scalar tensor theory of grav- 
itation |2l|. Other scalar field models include k-essence 
field [lj branes [22j and fluid models like the Chap- 
lygin gas model and its generalizations [23| . There are 
also some phenomenological models (25| , field theoretical 
and renormalization group based models (see e.g. [26j). 
models that unify dark matter and dark energy [28| , holo- 
graphic dark energy models [2t|; QCD dark energy [3(| 
and many others like those based on horizon thermody- 
namics (e.g. see HU). 
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Different models of dark energy which have the same 
background evolution are indistinguishable purely from 
distance measurements. Evolution of perturbations in 
these models is expected to break this degeneracy. The 
Integrated Sachs Wolfe (ISW) effect can distinguish a 
cosmological constant from other models of dark energy, 
especially ones with a dynamical dark energy jlq . Dark 
energy perturbations have been extensively studied in the 
linear approximation [U, 1321433 ] . Perturbations in dark 
energy affect the low I quadrupole in the CMB angular 
power spectrum through the ISW effect 0, HH ■ It was 
shown in Ref.[24| that dark energy perturbations affect 
the low I quadrupole in the CMB angular power spectrum 
through the ISW effect. For models with w > — 1 this 
effect is enhanced while for phantom like models it is sup- 
pressed. In these models dark matter perturbations and 
dark energy perturbations are anti-correlated for large 
effective sound speeds. This anti-correlation is a gauge 
dependent effect [HJ . There are several other studies of 
perturbations in dark energy (35|, including some that 
deal with evolution of spherical perturbations [3(1 l37l |. 
For canonical scalar field dark energy, the perturbations 
in matter are enhanced by the presence of dark energy 
perturbations in comparison with smooth dark energy 
model 40]. The matter perturbations in fluid models 
are suppressed compared to corresponding homogeneous 
dark energy scenario [4l| . As long as the speed of prop- 
agation of perturbations 'cj' is positive the evolution of 
matter perturbations is indistinguishable from a smooth 
dark energy model. This is true for scales smaller then 
the Hubble radius. Dark energy perturbations in a fluid 
model with an appropriate c 2 s emulate that of a scalar 
field model very well below the Hubble scale but start to 
differ at larger scales. Therefore the fluid model is not 
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a good approximation at these scales. This also implies 
that the growth of perturbations at large scales depends 
on the details of the model even though the background 
evolution is the same. A separate analysis is therefore 
required for every model. 

In this paper we consider different scalar held models 
to study evolution of dark energy perturbations. We con- 
sider two different types of potentials, classified as 'thaw- 
ing' and 'freezing'' in Ref. [39] . For potentials with a 
thawing behavior, the scalar field is frozen at early times 
and starts to roll down the potential at late time. Hence 
the equation of state of dark energy starts near w = — 1 
at early times to w > — 1 at late times. In contrast, in 
the case of potentials with freezing behavior, the scalar 
field rolls down the potential and approaches the mini- 
mum of the potential, with the equation of state going 
from w > — 1 to freezing at w = —1 at late times. Due 
to the different evolution of the equation of state, it is 
expected that the perturbations in dark energy evolve 
differently in these classes. For freezing potentials, at 
early times, the equation of state deviates from w = — 1 
hence the perturbations in dark energy are expected to 
be enhanced. Whereas, in thawing potential scenarios, 
dark energy perturbations become significant only at late 
times. Since dark energy perturbations enhance pertur- 
bations in nonrelativistic matter, models with freezing 
type potentials will have more enhanced perturbations 
in matter than those with thawing type behavior. This 
can be an additional tool (apart from distance measure- 
ments) to distinguish between these types of models. 

The paper is organized as follows. In Section [TT] we 
discuss evolution of perturbations in scalar field poten- 
tials in the thawing and freezing classes. The results are 
summarized in the concluding Section [ill! 



II. 



PERTURBATIONS IN SCALAR FIELD 
DARK ENERGY 



where dot denotes derivative with respect to the coordi- 
nate time t and vnr is the potential for the matter pecu- 
liar velocity such that 8v,i = ViVjvij. In these equations, 
we have Fourier decomposed the perturbed quantities 
such as 6<j>, 5pnr and vnr. and replaced V 2 by —k 2 , 
where k is the wave number defined as k — 2tt/X with 
A being the comoving length scale of perturbation. In 
these equations all the perturbed quantities correspond 
to the amplitude of perturbations in the k th mode. 

We assume a flat spatial cosmology and we choose the 
initial equation of state w to be very close to — 1 at early 
times (say redshift of Z{ = 100). The initial value of the 
scalar field needs to be fine tuned such that universe be- 
gins to accelerate at late times. We fine tune the param- 
eters such that the matter density parameter is within 
the range allowed by the present day observations. For a 
system of pressureless matter and a scalar field, the dy- 
namics of perturbations is uniquely determined by $(i) 
and 8<p(t). Hence we need two second order equations 
which connect $(i) and 8<j>{t). As one of these equations, 
we choose the third equation in system[21 The dynamical 
equation for the perturbations in the scalar field 8<f>(t) is 
given by : 
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■2<^V'{(j))-^ + V"{(t>)8(j) = Q (3) 



These equations complete the system and we can then 
calculate the fractional density perturbation for dark en- 
ergy and for nonrelativistic matter as 
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To describe dark energy perturbations, we choose the 
Newtonian gauge. We choose to work in this gauge as in 
this case we can directly relate the metric perturbation 
$ to the gravitational potential perturbation. In the ab- 
sence of anisotropic stress, the perturbed metric can be 
written in the form 

ds 2 = (1 + 2<$>)dt 2 - a 2 {t) [(1 - 2<5>)8 aP dx a dx p ] (1) 

where $ is the gauge invariant potential defined in [38| . 

The linearized Einstein equations obtained from the 
above metric are given by 
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To solve the equations numerically, we introduce the 
following dimensionless variables 
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Here $at is the normalized gravitational potential with 
$i n being the value of the metric potential at the initial 
time t — ti n . In terms of these two new variables, the 
equations are 



8 , 

4-$; 

s 

1 

2 



N 



N 



(7) 



y'Sy' - <P N y' 2 + V"(y)8y 



Sy" + 3-Sy'-W N y' 



k 2 Sy 



s s' 
+ 2$ N V'(y) + V"(y)y' ■■ 



(8) 



3 




V(0)=M 4 exp(-a 1/2 0) 
a = 0.1 




0.1 

a/ a 



o 



FIG. 1: The figure shows evolution of nonrelativistic matter density contrast as a function of the scale factor for the exponential 
potential (TH1). The thick solid line is the evolution of density contrast in the cosmological constant model. The solid line 
is for a homogeneous dark energy scenario whereas the dashed line show the evolution if dark energy perturbations are taken 
into account. 



Here k — kc/Ho, is the wave number scaled with respect 
to the Hubble radius, Hq being the Hubble parameter 
at the present epoch. The prime denotes derivative with 
respect to the variable tHo. 

We assume that the perturbation in the scalar field in 
the matter dominated epoch is negligibly small compared 
to other perturbed quantities such as $ and Snr- Hence 
we can treat the scalar field to be initially homogeneous. 
This corresponds to setting the initial condition 5y ln = 
and Sy' in — 0. The only initial condition that needs to be 
determined is the value of 3>jy. at t — t m . In the matter 

dominated epoch, the gravitational potential <I>(i) = 
for all values of the wave number fc, therefore we can set 
the initial condition &' Nia (k) = 0. 

If dark energy is a cosmological constant, then it does 
not cluster and the gravitational field $ decays when the 
cosmological constant dominates. For w ^ — 1 too, the 
potential $ remains constant in the matter dominated 
era and starts to decay when dark energy contribution 
becomes significant. This rate of decay is dependent on 
the value of w. For a canonical scalar field, dark en- 
ergy perturbations are correlated with the matter per- 
turbations and they enhance matter perturbations [4Cj . 
If dark energy is a barotropic fluid, the matter perturba- 
tions are suppressed as compared to the corresponding 
homogeneous model [4l|. This difference is due to the 
fact the homogeneous limit is achieved differently for the 



two models. For a scalar field, even if one ignores spa- 
tial gradients there still remains a contribution to SP. 
This contribution cancels with a corresponding contri- 
bution from the pressure term due to the background. 
This term cancels for homogeneous scalar field dark en- 
ergy and leads to a suppression in matter perturbations. 
For fluid models SP vanishes and the residual pressure 
term due to the background evolution makes the evolu- 
tion different from that in a scalar field. Therefore, as- 
suming dark energy to be homogeneous leads to a large 
difference in matter density contrast. 

For sub-Hubble scales, a fluid model effectively emu- 
lates a scalar field model [4l|. At scales smaller then the 
Hubble radius, dark energy can be assumed to be homo- 
geneous. It is sufficient to assume a fluid approximation 
and the number counts of clusters at different redshifts 
15] distinguish deviation from the cosmological constant 
36]. At larger scales, where dark energy perturbations 
may play a significant role, the fluid analogy breaks down 
and the evolution of matter density contrast depends on 
individual scalar field models. This also implies that the 
growth of perturbations at large scales depends on the de- 
tails of the model even though the background evolution 
is the same. A separate analysis is therefore required for 
every model. Since dark energy perturbations are signifi- 
cant at large scales, i.e., scales larger then Hubble radius, 
we will restrict the subsequent discussion to large scales. 
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FIG. 2: The plot on the left shows evolution of matter density contrast for the concave potentials (TH2) with respect to the 
scale factor. The figure on the left shows evolution of nonrelativistic matter density contrast as a function of the scale factor 
for the linear potential. The figure on the right shows evolution of nonrelativistic matter density contrast as a function of the 
scale factor for the quadratic potential. The line styles are the same as in Fig. 1. 



In particular, the reference scale we choose here is the 
length scale A = 10 5 Mpc. 

The scalar field models have been classified as freez- 
ing or thawing in [39[ . The classification is based on the 
difference in the way equation of state of dark energy 
evolves. The freezing type models are fast roll models 
with a steep potential. The field remains subdominant 
and at late times becomes dominant and drives the ac- 
celeration of the expansion of the universe. The thawing 
models have a potential which is nearly flat. Therefore, 
the equation of state starts with the value w = —1 at 
early times and the energy density of the field is negligi- 
ble. As the universe expands, the energy density of the 
field becomes comparable to the background energy den- 
sity. The equation of state of dark energy deviates away 
from its frozen value, hence the scalar field thaws. 

The dark energy equation of state evolves differently 
in these classes of models. Therefore, the evolution of 
perturbations in freezing type and thawing type models 
is expected to be different. Using present day distance 
measurement observations, it is not possible to distin- 
guish between different dark energy models from The 
ACDM model and also between various scalar field mod- 
els. These observations include Present day Supernova 
observations and data from the Baryon Acoustic Oscil- 
lations measurements [27| . This degeneracy is expected 
to be removed if dark energy perturbations are included. 



The difference in these models is expected to make a 
significant contribution to integrated effects such as the 
ISW effect [Hj]. 

To study scalar field perturbations, we consider the 
following potentials: 

1. For 'thawing' behavior 

(a) TH1: Exponential potential : 
V(4) = M 4 exp(-^(j)/Mp) [HI 

(b) TH2: Polynomial (concave) potential : 

v{<j>) = M 4 ~ n (j) n dnii 

2. For 'freezing' behavior 

(a) FR1: Inverse power potential : 
V((j)) = M 4 + n (f>- n di 

(b) FR2: V(4>) = M i+n <f>- n exp{a4> 2 / 'M%) pj] 

We have labeled the scalar field potentials as TH (thaw- 
ing) and FR (freezing), and henceforth will be using 
these labels to refer to the models. Since Einstein equa- 
tions imply homogeneous dark energy to be inconsistent, 
we will compare the evolution of perturbations in scalar 
field dark energy with the cosmological constant. We 
have however also considered the homogeneous limit in 
all these models. 

The growth of perturbations in scalar field model with 
an exponential potential (TH1) has been studied in [40| . 



5 




V(0) = M 4+n /0 r 



n = 2 



0.1 

a/ a 



o 



FIG. 3: The figure shows evolution of nonrelativistic matter density contrast as a function of the scale factor for the inverse 
potentials (FR1) for the parameter values mentioned. 



For this potential, if a = 1, the present day equation of 
state is w = —0.86 and the acceleration of the universe 
starts at redshift z = 0.81. For smaller a the models w is 
smaller, for instance if a = 0.1, the equation of state at 
present is w = —0.95. In Fig. Q]we plot the evolution of 
matter perturbations at scale A = 10 5 Mpc as a function 
of the scale factor a. For smaller a where the present 
equation of state is closer to that of a cosmological con- 
stant, there is no significant difference if we assume dark 
energy to be homogeneous or if we assume dark energy to 
be clustered. If a = 1 there is a significant enhancement 
of matter perturbations [40j • Compared to the cosmolog- 
ical constant model, matter perturbations in this model 
are suppressed, with approximately 6% deviation. 

In Fig. [5] we show the evolution of matter density con- 
trast for the concave potential (TH2) with n = 1 and 
n = 2. The present day of equation of state is w « —0.94 
if we choose n — 1 and reaches w « —0.3 for n = 2. If 
the potential is quadratic, the gravitational potential ini- 
tially decays, and in the future starts to oscillates. This is 
due to the fact that the scalar field reaches the minimum 
of the potential and the equation of state for dark energy 
then oscillates between w = — 1 and w = 1. These oscil- 
lations are evident at earlier times at very large scales, 
say at A = 10 5 Mpc. At large scales, there is a significant 
enhancement in matter perturbation as compared to the 
case when dark energy is assumed to be homogeneous. In 
contrast, the gravitational potential $ continues to decay 
if n = 1. For n = 1 the matter perturbations are sup- 



pressed as compared to the potential with n = 2 but they 
continue to grow and become larger when oscillations set 
in the case of the quadratic potential. Matter pertur- 
bations in linear potential follow those in cosmological 
constant model very closely. In quadratic scalar field po- 
tential case, these perturbations are enhanced compared 
to the ACDM model at early times and are suppressed 
at late times. At a ~ 0.4, the percentage enhancement 
as compared to the cosmological constant model is w 5 
%. 

For potential FR1, the scalar field begins to roll down 
the potential at early times and begins to freeze at late 
times. In this freezing type evolution, the equation of 
motion of dark energy starts away from w = — 1 and ap- 
proaches this value as the scalar field slows down. The 
epoch at which this freezing occurs is subject to fine tun- 
ing of parameters. If the scalar field is frozen before the 
present time, the duration of the matter dominated phase 
is then too small, i.e., scalar field approaches the kinetic 
energy dominated phase before reaching z = 1000. To 
achieve a sufficiently prolonged matter dominated phase, 
the scalar field should freeze in far future. The dark 
energy equation of state reaches a maximum of approxi- 
mately —0.7, and we tune the initial value of <p such that 
we achieve Qnr ~ 0.3 at present. In Fig. [3] we show 
the evolution of matter density contrast with scale fac- 
tor. There is an enhancement in matter perturbations 
with an increasing value of n. We have plotted the den- 
sity contrast for n = 1 and for n = 2. In this model, since 
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FIG. 4: The figure shows evolution of nonrelativistic matter density contrast as a function of the scale factor for the potential 
V{4>) = M 4+n cj>-"exp(a<p 2 /M 2 P ) (FR2). The plot of the left corresponds to a = 0.6 and n = 4 . The plot on the right shows 
the density contrast for a = 0.8 and n = 4. 



the value of equation of state deviates away from w = — 1 
at early times, there is a significant enhancement in dark 
energy perturbations and hence in matter perturbations. 
With n — 1, perturbations in nonrelativistic matter are 
comparable to those in cosmological constant model, and 
with n = 2 they are enhanced by approximately 7 % at 
a m 0.3 with respect to the cosmological constant model. 

For the potential FR2, the freezing behavior is achieved 
earlier without trading away a viable duration of the mat- 
ter dominated era. For a given n, the present value of 
the equation of state moves further away from w = — 1 
as a increases. Also for a constant a, the value of the 
equation of state increases with an increase in n. As in 
all other models, at very early time we assume w = — 1. 
For instance, if a = 0.6, for n = 1 reaches w ~ —0.9975 
at redshift s» 0.4 and then starts to approach the freezing 
behavior. If a = 0.8 and n = 4, the maximum value the 
equation of state reaches is w = —0.82. In Fig. 01 we 
show the evolution of matter density contrast at A = 10 5 
for two different parameter sets and also in comparison 
to that in ACDM model. 

In Fig. [5] and Fig. [HI we show evolution of matter per- 
turbations in different models considered here at different 
length scales. At small scales, say at A = 50 Mpc, the 
evolution is indistinguishable from that in a homogeneous 
model of dark energy. The evolution of perturbations in 
different in different scalar field models but the limit of 
a smooth dark energy works well at small length scales. 



At this scale, except for the inverse power law potential 
model, evolution in the scalar field models closely fol- 
lows that in the cosmological constant model. The fluid 
approximation too works well at these scales. Explicit 
dependence of dark energy perturbations is significant at 
large scales. 

At scales larger than the Hubble radius, the evolution 
of perturbed quantities are significantly different from 
those in homogenous dark energy model. In Fig. [51 we 
plot the gravitational potential and matter density con- 
trast at length scales A = 10 4 and A = I0 5 Mpc. The 
freezing models of scalar field, at early times have sig- 
nificant enhancement in matter perturbations at early 
times. There is a significant enhancement of matter per- 
turbations as compared to the other models. This en- 
hancement is more for models with higher n. Among the 
scalar field models considered here, models FR1 have the 
highest rate of growth of perturbations. The dashed line 
denotes FR2, with n — 2 and a = 0.8. As compared 
to the cosmological constant, FR2 has a higher rate of 
growth of nonrelativistic matter density contrast. The 
dot dashed line corresponds to the potentials TH2 with 
n = 2. However, the quadratic potential shows a larger 
growth rate at early times and then shows a downward 
trend when the scalar field approaches the minimum of 
the potential. At early times, TH2 has more enhance- 
ment in matter density contrast than the FR2 potential. 
This is because the dark energy equation of state devi- 




FIG. 5: This figure shows evolution of the gravitational potential $ and matter density contrast Snr as a function of scale 
factor for the four different models at length scale A = 50 Mpc. The solid line shows the evolution for the scalar field potential 
V{4>) = M i+n <t>~ n (FR1) with n = 2. The dashed line is for the potential V{4>) = M 4+n exp{a<t> 2 / 'Mj>) (FR2) with 
parameters a = 0.8 and n — 4. The dotted line corresponds to the thawing exponential potential V((f>) = M i exp{—y/a.(f>/Mp) 
(TH1) where a — 1 and the dot dashed line corresponds to the polynomial potential V((j>) = M 4 ~ n n (TH2) for n — 2. At 
early times, the freezing potentials have matter perturbations enhanced as compared to the thawing potentials. The evolution 
at this scale is same for a homogeneous dark energy and for a clustering dark energy. The thick solid line corresponds to the 
ACDM model. 



ates from that of a cosmological constant for the poly- 
nomial potential (TH2) more than for the inverse power 
exponential potential (FR1). At late times, for TH2 the 
scalar field reaches the minimum of the potential and the 
equation of state oscillates between —1 and +1 and the 
perturbations in dark energy and consequently matter 
perturbations begin to oscillate in future. The late time 
suppression shown in the figure is due to the beginning 
of the oscillatory phase. In the case of an exponential 
potential (TH1), matter density contrast is suppressed 
compared to all other models. 



III. CONCLUSIONS 

In this paper, we analyze the growth of perturbations 
in scalar held dark energy scenarios. The assumption 
that the distribution of dark energy (with w ^ —1) is 
homogeneous at all length scales is inconsistent with the 
equivalence principle. On length scales comparable to 
or greater than the Hubble radius, the perturbations in 
dark energy can become comparable to perturbation in 
matter if u>de ^ — 1- For scales smaller than the Hubble 
radius, perturbations in dark energy can be neglected in 



comparison with the perturbation in matter at least in 
the linear regime. Hence any deviations from the cosmo- 
logical constant model can be explored by assuming dark 
energy as a homogeneous component and the scalar field 
models can be approximated well by parameterized fluid 
models. However, on much larger scales, if the equation 
of state parameter deviates from -1, then perturbations 
in dark energy do influence matter power spectrum. 

For canonical scalar fields, a clustering dark energy 
enhances matter perturbations as compared to the cor- 
responding homogeneous dark energy scenario. This en- 
hancement more pronounced at scales larger than the 
Hubble radius. In particular, we study two broad classes 
of canonical scalar field models, namely the thawing and 
freezing models. The equation of state of dark energy 
evolves differently in these two classes of models. This af- 
fects the growth of perturbations in these different types. 
For fast rolling, freezing type models, the dark energy 
equation of state deviates away from that of a cosmolog- 
ical constant at early times and freezes to w = — 1 at 
late times. In slow roll thawing models, the scalar field 
remains at a constant w = — 1 and starts to deviate away 
from this value at late times. The present day observa- 
tions, which are based on distance measurements, cannot 
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FIG. 6: The upper panel in this figure shows evolution of the gravitational potential $ and matter density contrast Snr as 
a function of scale factor for the four different models at length scale A = 10 4 Mpc. The line styles used for different models 
are same as in Fig. 5. At early times, the freezing potentials have matter perturbations enhanced as compared to the thawing 
potentials. The lower part of the figure shows the evolution of gravitational potential and matter density contrast at length 
scale A = 10 5 . 
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distinguish between these models. 

We studied evolution of matter perturbations in the 
presence of a clustering dark energy. The models we have 
considered in this paper are within the range allowed by 
distance measurement observations. For corresponding 
homogeneous dark energy models and clustering dark en- 
ergy models, there are significant changes in the matter 
density contrast evolution. All the canonical scalar field 
models studied here show an enhancement in matter per- 
turbations if dark energy is perturbed. Although we have 
not considered phantom like models in this paper, it is 
worth mentioning that in these models, dark energy per- 
turbations suppress matter perturbations. 

In general, freezing type models have a higher rate of 
growth of density contrast at early times. This is due 
to the fact the in these models, the equation of state of 
dark energy is further away from that of a cosmological 
constant. Hence dark energy perturbations play a promi- 
nent role at earlier times. For thawing type models, in 
general, dark energy perturbations affect the matter per- 
turbations at late times. For most of the evolution, the 
matter perturbations remain close to those in cosmolog- 
ical constant model and being to deviate as the field be- 



gins to thaw. There are significant deviations in the way 
density contrast grows, not only between various mod- 
els but also from the concordant cosmological constant 
model. Apart from different way matter perturbations 
grow in the freezing and thawing classes of models, mod- 
els within the same class also have variation in the be- 
havior of the density contrast. Observable changes in the 
angular power spectrum at large scales are limited by the 
cosmic variance and therefore CMB data alone will be in- 
sufficient to distinguish between these models. Since the 
scales at which dark energy perturbations are relevant 
are large, the primary contribution to CMB anisotropics 
is through the ISW effect. Therefore it is important to 
study the ISW cross correlation with large scale structure 
indicators. 
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